Analysis of Tensile Strength and Failure Mechanism Based on Parallel Homogenization Model for Recycled Concrete

In this paper, a parallel homogenization model for recycled concrete was proposed. A new type of finite element method, the base force element method, based on the complementary energy principle and the parallel homogenization model, is used to conduct meso-level damage research on recycled concrete. The stress–strain softening curve and failure mechanism of the recycled concrete under uniaxial tensile load are analyzed using the nonlinear damage analysis program of the base force element method based on the parallel homogenization model. The tensile strength and destructive mechanisms of recycled concrete materials are studied using this parallel homogenization model. The calculation results are compared with the results of the experiments and meso-level random aggregate model analysis methods. The research results show that this parallel homogenization analysis method can be used to analyze the nonlinear damage analysis of recycled concrete materials. The tensile strength, stress–strain softening curve, and crack propagation process of recycled concrete materials can be obtained using the present method.


Introduction
Recycled aggregates are aggregates made from waste concrete through a series of processing methods. Recycled aggregate concrete is concrete made by replacing part, or all, of the natural aggregate with recycled aggregate. It has been widely valued as a green and environmentally friendly building material [1].
Many scholars have carried out a lot of experimental research on recycled concrete, and some research results have been obtained [1][2][3][4][5]. However, the test cycle is long, the cost is high, and it is difficult to measure the internal stress, strain, and failure mechanisms of the material. Therefore, it is very useful to carry out numerical simulation research on recycled concrete. In reference [6], a plastic-damage constitutive models are employed in numerical studies on recycled concrete under uniaxial compression and uniaxial tension loadings to predict the overall mechanical behavior, particularly the stress-strain relationship. In reference [7], a statistical analysis on its composition has been performed considering the randomness in properties of old adhered mortar around recycled aggregate. Peng et al. [8][9][10][11] proposed the base force element method and used this new type of finite element method to carry out a numerical simulation analysis on the recycled concrete, and conducted uniaxial tensile and compression tests under static and dynamic loads to study its mechanics and performance. In 2016, Rajendra [12] established a virtual crack model and a double-K fracture model, and determined the fracture parameters of recycled concrete with different coarse aggregate contents. In reference [13], a stochastic elastic FEM analyses model was established based on the Nano-indentation technique for recycled concrete at three different scales to obtain the ef-fective elastic moduli and Poisson's ratios, also the correlations of recycled concrete were studied. Anuruddha [14] investigated the influence of the old mortar content on the elastic module and the strength of recycled concrete, and found that the pressure strength of recycled concrete decreased with the increase in mortar content. Job et al. [15] used numerical simulation and regression analysis methods to study various mechanical properties of recycled concrete, and obtained the influence of different replacement rates of recycled aggregate on the strength of recycled concrete. In 2019, Tan [16] carried out a two-dimensional numerical simulation of recycled concrete based on the discrete element method, and mainly studied the influence of the weak link between the old and new interfaces on the damage and failure process of recycled concrete. Guo [17] has developed a creep coupling model for the heterogeneity of recycled concrete, to investigate the influence of recycled aggregate on the creeping of recycled concrete. Kazemian et al. [18] conducted experimental research on the compressive strength, flexural strength, and fracture energy of recycled concrete, and compared the mechanical properties of treated and untreated recycled concrete. There are also some scholars who have carried out research works in this field or other types of concrete [19][20][21][22][23][24][25][26].
In the paper, a homogenization analysis method will be used to establish a parallel homogenization model for recycled concrete materials. A new type of finite element method, the base force element method based on the complementary energy principle, is used to conduct meso-level damage research on recycled concrete. The stress-strain softening curve and failure mechanisms of recycled concrete under uniaxial tensile load are analyzed using the nonlinear damage analysis program of the base force element method.

Random Aggregate Model
The main difference between natural concrete and recycled concrete is that the outer layer of recycled concrete aggregate is attached with a layer of old mortar. The random aggregate model represents the recycled concrete in the form of each phase medium, so as to facilitate the subsequent mesh division and calculation at the meso level. The circular aggregate model was adopted to simplify the recycled aggregate into two concentric circles, as shown in Figure 1. The macro-mechanical properties of the whole structure are obtained through the analysis and calculation of the random aggregate model. It is the link between the macro-structure and the macro-mechanical properties of recycled concrete. to analyze the stress-strain relationship of recycled concrete under different replacement rates. Anuruddha [12] investigated the influence of the old mortar content on the elastic module and the strength of recycled concrete, and found that the pressure strength of recycled concrete decreased with the increase in mortar content. Job et al. [13] used numerical simulation and regression analysis methods to study various mechanical properties of recycled concrete, and obtained the influence of different replacement rates of recycled aggregate on the strength of recycled concrete. In 2019, Tan [14] carried out a twodimensional numerical simulation of recycled concrete based on the discrete element method, and mainly studied the influence of the weak link between the old and new interfaces on the damage and failure process of recycled concrete. Guo [15] has developed a creep coupling model for the heterogeneity of recycled concrete, to investigate the influence of recycled aggregate on the creeping of recycled concrete. F. Kazemian et al. [16] conducted experimental research on the compressive strength, flexural strength, and fracture energy of recycled concrete, and compared the mechanical properties of treated and untreated recycled concrete. There are also some scholars in this field who have carried out research work on other types of concrete [17][18][19][20][21][22][23][24].
In the paper, a homogenization analysis method will be used to establish a parallel homogenization model for recycled concrete materials. A new type of finite element method, the base force element method based on the complementary energy principle, is used to conduct meso-level damage research on recycled concrete. The stress-strain softening curve and failure mechanisms of recycled concrete under uniaxial tensile load are analyzed using the nonlinear damage analysis program of the base force element method.

Random Aggregate Model
The main difference between natural concrete and recycled concrete is that the outer layer of recycled concrete aggregate is attached with a layer of old mortar. The random aggregate model represents the recycled concrete in the form of each phase medium, so as to facilitate the subsequent mesh division and calculation at the meso level. The circular aggregate model was adopted to simplify the recycled aggregate into two concentric circles, as shown in Figure 1. The macro-mechanical properties of the whole structure are obtained through the analysis and calculation of the random aggregate model. It is the link between the macro-structure and the macro-mechanical properties of recycled concrete. In Figure 1, the recycled concrete is treated as a 5-phase medium, including the aggregate, the old mortar, the new mortar, the old interface and the new interface. The aggregate center position is automatically generated by a computer program according to the Monte Carlo method. The aggregate size and particle number are calculated according In Figure 1, the recycled concrete is treated as a 5-phase medium, including the aggregate, the old mortar, the new mortar, the old interface and the new interface. The aggregate center position is automatically generated by a computer program according to the Monte Carlo method. The aggregate size and particle number are calculated according to the grading of recycled concrete. The placement of aggregate should ensure that each aggregate cannot overlap. The calculation and analysis are carried out on the basis of the random aggregate model of recycled concrete. First, the two-dimensional random aggregate is generated. After that, the mesh is divided and mapped to the model, the position of the element node is judged, and the attributes of each element are assigned. Aggregates are divided into coarse aggregates and fine aggregates. Fine aggregates refer to aggregate particles with a particle size of less than 5 mm. The influence of the fine aggregate is ignored, and the fine aggregate and mortar are regarded as a whole in the research of this article. The data of the particle size range of the coarse aggregate are obtained based on the experimental data.

Parallel Homogenization Model
A parallel method was used to homogenize the heterogeneous elements of recycled concrete, based on the basic idea of equivalent model of meso-damage element, and the equivalent model of meso-damage was established. The validation and numerical calculation of homogenized equivalent model are based on random aggregate model. A simple material strength equivalent method was adopted using the Kelvin-Voigt parallel model, ignoring transverse deformation, as shown in Figure 2. σ is the stress and E is the modulus of elasticity.
aggregate cannot overlap.
The calculation and analysis are carried out on the basis of the random model of recycled concrete. First, the two-dimensional random aggregate is After that, the mesh is divided and mapped to the model, the position of the el is judged, and the attributes of each element are assigned. Aggregates are d coarse aggregates and fine aggregates. Fine aggregates refer to aggregate part particle size of less than 5 mm. The influence of the fine aggregate is ignored, a aggregate and mortar are regarded as a whole in the research of this article. the particle size range of the coarse aggregate are obtained based on the ex data.

Parallel Homogenization Model
A parallel method was used to homogenize the heterogeneous elements concrete, based on the basic idea of equivalent model of meso-damage eleme equivalent model of meso-damage was established. The validation and nume lation of homogenized equivalent model are based on random aggregate mode material strength equivalent method was adopted to solve the Kelvin-Vo model, ignoring transverse deformation, as shown in Figure 2. σ is the stre the modulus of elasticity. The equilibrium equation of element stress is as follows: where A is the area of an equivalent element, A1 and A2 are the areas of tw media elements, respectively. The strain compatibility equation is as follows: The mean stress-mean strain relationship is as follows: It can be obtained from the following formulas: The equilibrium equation of element stress is as follows: where A is the area of an equivalent element, A 1 and A 2 are the areas of two different media elements, respectively. The strain compatibility equation is as follows: The mean stress-mean strain relationship is as follows: It can be obtained from the following formulas: For the composite recycled concrete model, it can used the volume fraction to calculate , then we can obtain the following formula: where c 1 and c 2 are the percentage of the two different media elements to the total volume of the equivalent elements, respectively. The equivalent Poisson's ratio is as follows:

Finite Element Meshing and Homogenization
The finite element method is the process of dividing a continuous object into simple and regular elements. It is an effective method for calculating and analyzing objects with heterogeneous materials. Scholars have proposed many meshing methods, as the finite element method is widely used in the scientific community. The most common elements are as follows: one-dimensional rod element, two-dimensional triangle element, quadrilateral element and three-dimensional tetrahedron element, pentahedron element and hexahedron element. The finer and denser the mesh, the more accurate the calculation result for the traditional finite element method. However, at the same time, the more computational elements, the lower the computational efficiency.
In this paper, the meshing method of the two-dimensional quadrilateral element is selected based on the base force element method of complementary energy principle. The mesh is divided according to the size of the specimen and the element size. The midpoint of each side of the quadrilateral element is used as the calculation point, as shown in Figure 3. Then the mesh is mapped to the random aggregate model.
( 1 2 1 c c + = ), then we can obtain the following formula: where c1 and c2 are the percentage of the two different media elements to the total of the equivalent elements, respectively.
The equivalent Poisson's ratio is as follows:

. Finite Element Meshing and Homogenization
The finite element method is the process of dividing a continuous object into and regular elements. It is an effective method for calculating and analyzing obje heterogeneous materials. Scholars have proposed many meshing methods, as t element method is widely used in the scientific community. The most common e are as follows: one-dimensional rod element, two-dimensional triangle element, lateral element and three-dimensional tetrahedron element, pentahedron elem hexahedron element. The finer and denser the mesh, the more accurate the cal result for the traditional finite element method. However, at the same time, the mo putational elements, the lower the computational efficiency.
In this paper, the meshing method of the two-dimensional quadrilateral ele selected based on the base force element method of complementary energy princi mesh is divided according to the size of the specimen and the element size. The m of each side of the quadrilateral element is used as the calculation point, as show ure 3. Then the mesh is mapped to the random aggregate model. In Figure 4, the fine mesh of the random aggregate meso-model is divided f then the coarse mesh is formed by homogenization. The division method and ju rule of coarse mesh are the same as those of fine mesh, and the judgment rule mined according to the location of element nodes. The attribute of the element is as aggregate (or old mortar, or new mortar), when more than or equal to 3 of nodes of an element are projected on aggregate (or old mortar, or new mortar) m The element is defined as old interface element, when some element nodes fall i gate medium and some fall in old mortar. Similarly, the element is defined as interface when some nodes of an element fall on the old mortar medium and som fall on the new mortar medium. After determining the attributes of fine eleme number of small elements of each attribute was counted, the area was calculated ing to the element size, and the proportion of each component in the coarse mesh In Figure 4, the fine mesh of the random aggregate meso-model is divided first, and then the coarse mesh is formed by homogenization. The division method and judgment rule of coarse mesh are the same as those of fine mesh, and the judgment rule is determined according to the location of element nodes. The attribute of the element is judged as aggregate (or old mortar, or new mortar), when more than or equal to three of the four nodes of an element are projected on aggregate (or old mortar, or new mortar) medium. The element is defined as old interface element, when some element nodes fall in aggregate medium and some fall in old mortar. Similarly, the element is defined as the new interface when some nodes of an element fall on the old mortar medium and some nodes fall on the new mortar medium. After determining the attributes of fine elements, the number of small elements of each attribute was counted, the area was calculated according to the element size, and the proportion of each component in the coarse mesh element was calculated. The coarse mesh is equivalent to uniform single-attribute element by homogenization method.
Materials 2022, 15, x FOR PEER REVIEW was calculated. The coarse mesh is equivalent to uniform single-attribute elem mogenization method. The mesh division and attribute assignment of the elements were prog calculate, and the mesh node numbers and coordinates of all the specimens we using Fortran language program. The attribute distribution of an element c tained from the element attribute file by the calculating program. The compo heterogeneous element containing heterogeneous medium could be obtaine are provided for the subsequent homogenization calculation.
The divided elements are mapped to the random aggregate model and large element is subdivided to obtain small elements. The proportions of the a the small elements, the number of components, and equivalent parameter calc determined. The equivalent parameters are assigned to large elements.
The distribution of equivalent elastic modulus in the model is obtained. T lent elastic modulus distribution of the parallel equivalent method of the test obtained, as shown in Figure 5. In the red area of Figure 5, all elements are aggregate. Therefore, the com tomatic judgment belongs to the same medium. Computers do not use equ cessing, that is, the elastic modulus is the elastic modulus of aggregate. The mesh division and attribute assignment of the elements were programmed to calculate, and the mesh node numbers and coordinates of all the specimens were obtained using Fortran language program. The attribute distribution of an element could be obtained from the element attribute file by the calculating program. The composition of the heterogeneous element containing heterogeneous medium could be obtained. The data are provided for the subsequent homogenization calculation.
The divided elements are mapped to the random aggregate model and then each large element is subdivided to obtain small elements. The proportions of the attributes of the small elements, the number of components, and equivalent parameter calculations are determined. The equivalent parameters are assigned to large elements.
The distribution of equivalent elastic modulus in the model is obtained. The equivalent elastic modulus distribution of the parallel equivalent method of the test specimen is obtained, as shown in Figure 5.  The mesh division and attribute assignment of the elements were programmed to calculate, and the mesh node numbers and coordinates of all the specimens were obtained using Fortran language program. The attribute distribution of an element could be obtained from the element attribute file by the calculating program. The composition of the heterogeneous element containing heterogeneous medium could be obtained. The data are provided for the subsequent homogenization calculation.
The divided elements are mapped to the random aggregate model and then each large element is subdivided to obtain small elements. The proportions of the attributes of the small elements, the number of components, and equivalent parameter calculations are determined. The equivalent parameters are assigned to large elements.
The distribution of equivalent elastic modulus in the model is obtained. The equivalent elastic modulus distribution of the parallel equivalent method of the test specimen is obtained, as shown in Figure 5. In the red area of Figure 5, all elements are aggregate. Therefore, the computer's automatic judgment belongs to the same medium. Computers do not use equivalent processing, that is, the elastic modulus is the elastic modulus of aggregate.
In the blue area in Figure 5, all elements are new mortar. Therefore, after the automatic judgment of the computer, the area belongs to the same medium. Computers do not perform the equivalent, that is, the elastic modulus is the elastic modulus of new mortar.
If all small elements inside the large element (equivalent elements) are old mortar, the area belongs to the same medium. There is no equivalent processing, that is, the elastic modulus is the elastic modulus of old mortar. This is shown in purple in Figure 5. In the red area of Figure 5, all elements are aggregate. Therefore, the computer's automatic judgment belongs to the same medium. Computers do not use equivalent processing, that is, the elastic modulus is the elastic modulus of aggregate.
In the blue area in Figure 5, all elements are new mortar. Therefore, after the automatic judgment of the computer, the area belongs to the same medium. Computers do not perform the equivalent, that is, the elastic modulus is the elastic modulus of new mortar.
If all small elements inside the large element (equivalent elements) are old mortar, the area belongs to the same medium. There is no equivalent processing, that is, the elastic modulus is the elastic modulus of old mortar. This is shown in purple in Figure 5.
Equivalent treatment is required if small elements within large elements have aggregate elements and interface elements, or if large elements contain aggregate elements, mortar elements and interface elements. The old interface transition zone and new interface transition zone belong to equivalent elements, as shown in Figure 5. The elastic modulus varies in this region. There are different colors.

Damage Model
The stress-strain relationship of material is a very important and complex problem in the case of material damage. The constitutive nature of material is an old and still open question, starting with the papers of Hudson et al., . Even more recently, Ferretti [30] has shown that the meso-scale curves (mean stress/mean strain and damage curves) are not constitutive, while it is possible to identify constitutive laws at the micro scale.
In this paper, the multi-line stress-strain relationship is adopted for calculation because the stress tends to be highly nonlinear when approaching the peak value under uniaxial tensile load, due to the non-uniformity of the recycled concrete. The material is damaged due to stretching. By introducing the scalar damage variable D, the relationship between the effective strain of the damaged material and Cauchy stress is as follows: The elastic modulus after damage can be expressed by the initial elastic modulus, if the effect of damage on Poisson's ratio is neglected, as follows: where E represents the elastic modulus after damage, and E 0 represents the initial elastic modulus. Therefore, the damage elastic modulus of the five-phase medium in recycled concrete can be expressed as (aggregate (ag), old mortar (om), new mortar (m), old interface (oitz), new interface (itz)).
where ω is the residual compressive strength coefficient; ε 0 is the peak strain; η is the residual strain coefficient; λ is the elastic strain coefficient; δ is the elastic compressive strength coefficient; ξ is the limiting strain coefficient; µ is the residual tensile strength coefficient. Subscripts t represent the tension of the element. The volume fraction of each phase medium can be simplified into area fraction for the two-dimensional random aggregate model. Assuming that the size of the large mesh is a and the size of the small mesh is b(a > b), and there are n small mesh element attributes in the large mesh determined as aggregate, then the area fraction of aggregate is c 1 = nb 2 /a 2 .
Similarly, c 0 , c 1 , c 2 , c 3 and c 4 are used to represent the area fraction of new mortar, aggregate, old interface, old mortar and new surface, respectively.
Based on the strain compatibility equation, we obtain the following: The elastic modulus of the parallel equivalent element is as follows: The thickness of the new and old interface is small for recycled concrete. The mesh element is larger when the homogenized equivalent model is used. A single element may contain multiphase media when the element size is larger than the thickness of the old mortar. Therefore, this article will use recycled concrete as an equivalent to the following three phases: mortar (m), aggregate (ag) and equivalent element (em). The mesoequivalent model of stress-strain relationship for three-phase medium is established, where The formula of elastic damage model of parallel equivalent element of recycled concrete is as follows: Formula (15) is the calculation formula of equivalent elastic modulus at each stage. This formula can be obtained by combining the stages in the constitutive model.
Below, we will deduce the calculation formula of equivalent tensile strength of recycled concrete according to energy equivalence. See Formula (16)-(18) for details. From this formula, the equivalent tensile strength of recycled concrete can be calculated. The equivalent element is a homogeneous element. The stored total strain energy W is W = W ag + W em + W om (the sum of the equivalent element, new mortar, and aggregate), when the equivalent element reaches the equivalent tensile strength f t eq . Because the strain energy is W = 1 2 f eq t εdV, where ε is the element strain, the following is obtained: Substituting ε = f t E into the above formula, we obtain the following: The mean stress/mean strain and damage curves of equivalent element can be obtained after equivalence of different media, as shown in Figure 6.
The mean stress/mean strain and damage curves of equivalent element can b tained after equivalence of different media, as shown in Figure 6. In Figure 6, the stress-strain curve of homogenized equivalent material is a m broken line form. The y-coordinate is the equivalent strength and the x-coordinate i equivalent strain. After the peak point, there is a period of strength decline, mater damaged and destroyed..
The slope of the curve can be calculated according to Formula (14) and (15). Zhu et al. [29] assumes that the properties of each component material confor the Weibull distribution, and considers the non-linear characteristics of material inh geneity, thus proposing a random mechanical model. The density function of the We distribution is as follows: In Figure 6, the stress-strain curve of homogenized equivalent material is a multibroken line form. The y-coordinate is the equivalent strength and the x-coordinate is the equivalent strain. After the peak point, there is a period of strength decline, material is damaged and destroyed..
The slope of the curve can be calculated according to Formula (14) and (15). Zhu et al. [31] assumes that the properties of each component material conform to the Weibull distribution, and considers the non-linear characteristics of material inhomogeneity, thus proposing a random mechanical model. The density function of the Weibull distribution is as follows: where m determines the shape of the Weibull distribution density function, and it represents the uniformity of the medium; u represents a random variable satisfying the Weibull distribution. This paper considers the random distribution of materials in each phase of recycled concrete, and its material parameters obey the Weibull distribution. The value of each parameter is shown in Table 1.

Uniaxial Tensile Loading Model
The loading model of a cubic recycled concrete specimen is established by using the calculation program of the base force element method of the complementary energy principle in MATLAB programming. The specimen is selected as 100 mm × 100 mm × 100 mm to carry out the uniaxial tensile numerical simulation test. First, the cube model is simplified into a two-dimensional model with a cross-section size of 100 mm × 100 mm, and the loading model is shown in Figure 7. Vertical loading is adopted during loading. A static displacement-controlled loading condition is adopted step-by-step, with a displacement of 0.01 mm for each stage.

Uniaxial Tensile Loading Model
The loading model of a cubic recycled concrete specimen is established by using the calculation program of the base force element method of the complementary energy principle in MATLAB programming. The specimen is selected as 100 mm × 100 mm × 100 mm to carry out the uniaxial tensile numerical simulation test. First, the cube model is simplified into a two-dimensional model with a cross-section size of 100 mm × 100 mm, and the loading model is shown in Figure 7. Vertical loading is adopted during loading. A static displacement-controlled loading condition is adopted step-by-step, with a displacement of 0.01 mm for each stage.

Numerical Simulation Results
Different random numbers are selected to be put in the aggregate, in order to obtain specimens with different aggregate distribution positions in the random aggregate model. Three two-dimensional numerical models, with three groups of different aggregate distributions and the same aggregate particle number, are selected, as shown in Figure 8.

Numerical Simulation Results
Different random numbers are selected to be put in the aggregate, in order to obtain specimens with different aggregate distribution positions in the random aggregate model. Three two-dimensional numerical models, with three groups of different aggregate distributions and the same aggregate particle number, are selected, as shown in Figure 8. The basic force element method based on the principle of complementary energy is used to calculate and analyze the medium damage by the calculation program. The parallel equivalent homogenization model of the generated three random specimens is numerically simulated in a uniaxial tensile test. The calculation results of the three specimens, random aggregate results, and test data [30] are listed in Table 2. Meanwhile, the full stress-strain curve is drawn with strain as the abscissa and stress as the ordinate. The calculation result of the parallel equivalent is shown in Figure 9.
The recycled concrete specimen is in the elastic stage at the initial stage, as shown in Figure 9. The stress begins to grow slowly, and reaches about 80% of the specimen's ultimate strength. When the strain of the local element is greater than the residual strain, it begins to enter the state of damage. As the strain increases, the stress decreases until it reaches zero.  The basic force element method based on the principle of complementary energy is used to calculate and analyze the medium damage by the calculation program. The parallel equivalent homogenization model of the generated three random specimens is numerically simulated in a uniaxial tensile test. The calculation results of the three specimens, random aggregate results, and test data [32] are listed in Table 2. Meanwhile, the full stress-strain curve is drawn with strain as the abscissa and stress as the ordinate. The calculation result of the parallel equivalent is shown in Figure 9. The recycled concrete specimen is in the elastic stage at the initial stage, as show Figure 9. The stress begins to grow slowly, and reaches about 80% of the specimen's mate strength. When the strain of the local element is greater than the residual stra begins to enter the state of damage. As the strain increases, the stress decreases un reaches zero. ) specimen 1 specimen 2 specimen 3 Meso Random Aggregate Test Data [32] Figure 9. Uniaxial tensile stress-strain curve of parallel model. The QuickWin module in Fortran is adopted to display the different stages of element with different colors, in order to obtain the damage diagram of the num simulation calculation of recycled concrete clearly and simply. The equivalent ele will be set as the same yellow-green color, the mortar as orange, the aggregate as and the failure element is represented by the black block. The failure mode of the nu ical simulation specimen is observed, as shown in Figure 10.
The law of tensile failure of recycled concrete can be observed from the damage diagram of specimens in the process of uniaxial tensile loading. The local failure o first, when the element begins to be loaded under a certain strain. With the increa loading strain, cracks gradually spread through the whole specimen, accompanied b final failure of the specimen. It can be found that the location at which the cracks dev is generally the location where the aggregate is more concentrated in the failure diag The main reason for this is that the surrounding strength of the reclaimed aggreg lower, and it is easier to reach the destruction stage first. The equivalent element con a multiphase medium, and contains the old interface and the new interface with low sile strength. Therefore, the element with the cracks that appear first is the equivale ement part. The crack development direction of the tensile failure of the homogen specimen is perpendicular to the loading direction, and the failure state is basically a izontal crack, which is consistent with the failure state of the random aggregate m and is combined with the actual law. The recycled concrete specimen is in the elastic stage at the initial stage, as shown in Figure 9. The stress begins to grow slowly, and reaches about 80% of the specimen's ultimate strength. When the strain of the local element is greater than the residual strain, it begins to enter the state of damage. As the strain increases, the stress decreases until it reaches zero.
The QuickWin module in Fortran is adopted to display the different stages of each element with different colors, in order to obtain the damage diagram of the numerical simulation calculation of recycled concrete clearly and simply. The equivalent element will be set as the same yellow-green color, the mortar as orange, the aggregate as blue, and the failure element is represented by the black block. The failure mode of the numerical simulation specimen is observed, as shown in Figure 10. first, when the element begins to be loaded under a certain strain. With the increase in loading strain, cracks gradually spread through the whole specimen, accompanied by the final failure of the specimen. It can be found that the location at which the cracks develop is generally the location where the aggregate is more concentrated in the failure diagram.
The main reason for this is that the surrounding strength of the reclaimed aggregate is lower, and it is easier to reach the destruction stage first. The equivalent element contains a multiphase medium, and contains the old interface and the new interface with low tensile strength. Therefore, the element with the cracks that appear first is the equivalent element part. The crack development direction of the tensile failure of the homogenized specimen is perpendicular to the loading direction, and the failure state is basically a horizontal crack, which is consistent with the failure state of the random aggregate model, and is combined with the actual law.

Discussion
Using the homogenization model of recycled concrete proposed in this paper, the equivalent element is used to replace the tiny random aggregate model element, and the number of available elements is greatly reduced.
In this way, the computing speed is increased and the computer memory is reduced. Table 3 shows the comparison data. The law of tensile failure of recycled concrete can be observed from the damage state diagram of specimens in the process of uniaxial tensile loading. The local failure occurs first, when the element begins to be loaded under a certain strain. With the increase in loading strain, cracks gradually spread through the whole specimen, accompanied by the final failure of the specimen. It can be found that the location at which the cracks develop is generally the location where the aggregate is more concentrated in the failure diagram. The main reason for this is that the surrounding strength of the reclaimed aggregate is lower, and it is easier to reach the destruction stage first. The equivalent element contains a multiphase medium, and contains the old interface and the new interface with low tensile strength. Therefore, the element with the cracks that appear first is the equivalent element part. The crack development direction of the tensile failure of the homogenized specimen is perpendicular to the loading direction, and the failure state is basically a horizontal crack, which is consistent with the failure state of the random aggregate model, and is combined with the actual law.

Discussion
Using the homogenization model of recycled concrete proposed in this paper, the equivalent element is used to replace the tiny random aggregate model element, and the number of available elements is greatly reduced.
In this way, the computing speed is increased and the computer memory is reduced. Table 3 shows the comparison data. The homogenization model enlarges the size of the grid element of the calculation model, and the equivalent parameters of the homogenization element are obtained by using the parallel equivalent formula.
For the two-dimensional random aggregate recycled concrete model with a size of 100 mm × 100 mm, when the element size of the random aggregate model is 0.5 mm and the element mesh size of the homogenized model is 2 mm, the element mesh size increases by 4 times, the number of elements becomes 1/16 of the original, and the calculation time is reduced by about 300 times.
Obviously, the homogenized equivalent model can greatly save calculation time and improve calculation efficiency, which provides a new way for future numerical simulation analysis and calculation.

Conclusions
(1) The parallel equivalent stress-strain relationship of the homogenization model is derived. The multi-line damage model of recycled concrete materials is established by using the homogenization analysis method.
(2) The non-linear basic force element analysis software and the homogenization preprocessing software for the homogenization analysis of recycled concrete have been developed, based on the basic force element method of the complementary energy principle.
(3) A parallel-equivalent homogenization model was used to perform a numerical calculation and analysis on the uniaxial tensile test of recycled concrete. The stress-strain softening curve, and the damage and failure process were obtained.
(4) The feasibility and rationality of the model establishment are verified by comparing the results of this method with the experimental data. (5) The calculation efficiency of the homogenization model has greatly improved. The calculation efficiency of this method is much higher than that of the mesoscopic damage analysis method based on the random aggregate model, and it can guarantee a certain calculation accuracy. (6) The research work in this paper shows that the base force element method based on the complementary energy principle and the parallel homogenization model can be used to analyze the meso-structure and mechanical properties of recycled concrete. It has the characteristics of high computational efficiency and can be used as an effective meso-analysis method for recycled concrete. (7) In the future, we will study the base force element method based on the complementary energy principle for dynamic damage analysis and three-dimensional analysis of recycled concrete.